CRYSTAL ENERGY FUNCTIONS VIA THE CHARGE 
IN TYPES A AND C 

CRISTIAN LENART AND ANNE SCHILLING 

^N) Abstract. The Ram- Yip formula for Macdonald polynomials (at i = 0) provides a statistic 

T-H which we call charge. In types A and C it can be defined on tensor products of Kashiwara- 

^^ Nakashima single column crystals. In this paper we prove that the charge is equal to the 

^N (negative of the) energy function on affine crystals. The algorithm for computing charge is 

rH much simpler and can be more efficiently computed than the recursive definition of energy in 

C^ terms of the combinatorial i?-matrix. 

o 
m 

I— I 1. Introduction 

o 

rj The energy function of affine crystals is an important grading used in one-dimensional con- 

figuration sums [lOl [11] and generalized Kostka polynomials [131 HSJ HB]. It is defined by the 
action of the affine Kashiwara crystal operators through a local combinatorial rule and the 
i?-matrix. 






> 



From a computational perspective, the definition of the energy is not very efficient, as it 

involves a recursive definition of a local energy, and also the combinatorial i2-matrix, for which 

^ not in all cases efficient algorithms exist. This leads us to the role of the charge statistic, 

^» which can be calculated very efficiently, as it only involves the detection of descents and the 

computation of arm lengths of cells in Young diagrams. 

■^ Charge was originally defined in type A by Lascoux and Schiitzenberger [20j as a statistic 

t~^ on words with partition content. It is calculated by counting certain cycles in the given word; 

O see Section [3j Lascoux and Schiitzenberger showed that the charge can also be defined as 

, the grading in the cyclage graph, and used it to express combinatorially the Kostka-Foulkes 

r* polynomials, or Lusztig's (/-analogue of weight multiplicities [34j, based on their Morris recur- 

. ^H rence. In type A^ Nakayashiki and Yamada [37J analyzed the subtle combinatorial relationship 

/\ between charge and the i?-matrix, showing that the energy coincides with the charge. In [19] 

it was observed that the cyclage is related to the action of the crystal operator /o on a tensor 
product of type A columns (by the Kyoto path model, the latter can be identified with an affine 
Demazure crystal). Thus, the results of Lascoux-Schiitzenberger and Nakayashiki- Yamada are 
rederived in a more conceptual way. See also the work of Shimozono \45\ [l6] for a more extensive 
discussion of the combinatorics involved in [19j, in the more general context of tensor products 
of type A Kirillov-Reshetikhin (KR) crystals of arbitrary rectangular shapes, as opposed to 
only column shapes. Charge for KR crystals of rectangular shape (or Littlewood-Richardson 
tableaux) was also defined in [43j using cyclage. 



^ 
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Lecouvey |24| [25] extended two approaches to the Lascoux-Schiitzenberger charge, namely 
cyclage and catabohsm, to types B, C, and D. He thus defined a charge statistic on the 
corresponding Kashiwara-Nakashima (KN) tableaux |17) . But he was only able to relate his 
charge to the corresponding Kostka-Foulkes polynomials in very special cases, as the original 
idea of Lascoux-Schiitzenberger based on the Morris recurrence, which he pursued, has limited 
applicability in this case. 

In this paper we use a charge statistic coming from the Ram- Yip formula [30] for Macdonald 
polynomials P^{x; q, t) of arbitrary type j35j at t = 0. The terms in this formula correspond to 
certain chains of Weyl group elements which come from the alcove walk model (this was defined 
in [3l|32l|33], and was then developed in subsequent papers, including |lQ]). The statistic is 
defined on the mentioned chains, and describes the powers of q. In [30] it is shown that, in 
types A and C, the chains are in bijection with elements in a tensor product of KR crystals of 
the form B '^. It is also shown that, under this bijection, the above statistic can be translated 
into a statistic on the elements of the mentioned crystal, which we call charge. Thus, we have 

(1.1) Pf,{x;q,0)= Yl gcharge{6)^wt(6)_ 

6e-B''i'^®B^2'i(g)... 

In type A, one can rewrite this formula as an expansion of the Macdonald P-polynomials in 
terms of Schur polynomials sx{x) 

(1-2) P^{x;q,0)=Y,Ky^,'{q)sxix); 

A 

here Kxif^'{q) is the Kostka-Foulkes polynomial and A' denotes the transpose of the partition A. 



A generalization of (1.2) to simply-laced types was given in [12]; in types A and D, this result is 
sharpened in [42" Section 9.2] by replacing the Kostka-Foulkes polynomials with the correspond- 
ing one-dimensional configuration sums (which are generating functions for the energy). Both 



(1.1) in type A and (1.2) are expressed combinatorially in terms of the Lascoux-Schiitzenberger 



charge, whereas the type C charge given by (1.1) is a new statistic. It is worth noting that 



the main ingredient in these charge constructions is the quantum Bruhat graph [2], which first 
arose in connection to Chevalley multiplication formulas for the quantum cohomology of flag 
varieties. 

Due to Schur-Weyl duality, weight multiplicities and tensor product multiplicities are both 
given by the Kostka matrix in type A. In fact, the Kostka-Foulkes polynomials are equal to 
Lusztig's g-analogue K of weight multiplicities as well as the one-dimensional configurations 
sums X, which are g-analogues of tensor product multiplicities. For other types, such corre- 
spondences are only known for large rank (or in the stable case); see for example [2b\ 1^5] [SU] . 



As is evident from Equation (1.1) the charge we are concerned with in this paper is naturally 
related to g-analogues of tensor product multiplicities or the one-dimensional sums X. 

The goal of this paper is to show in an efficient, conceptual way that the charge in [30] 
coincides with the energy function on the corresponding tensor products of KR crystals. We 
focus on types A and C, and expect to extend these results to types B and D. There are 
several reasons for not yet addressing the case of arbitrary types. First of all, the elements 
of the crystals B '^ of classical types are represented in a very explicit way, by KN columns. 
In fact, concrete descriptions of certain KR crystals of exceptional type have only been given 
in special cases; see for example J13t IT6l [361 151) . Furthermore, the stable one-dimensional 
configuration sums studied in [28] are all of classical type. However, we do expect that the 
main result of this paper would generalize to arbitrary types, if we were to phrase it in the 
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context of the alcove walk model and the statistic in the Ram- Yip formula mentioned above 



(see Section 7.3) 



We use the recent reinterpretation in [12] of the (global) energy function as the affine grading 



on a tensor product of KR crystals under Demazure arrows (see Definition 2.7). In type A, KR 
crystals are perfect and hence, by the Kyoto path model |14 | ll5 j . can be realized as Demazure 
crystals. By the result of [l2], Demazure arrows change energy by 1. Together with the result 
that charge is well-behaved under crystal operators, this proves the equality of energy and 
charge. For type C, we use the same approach, but in this case KR crystals are not perfect. 
There is still an embedding of a tensor product of KR crystals into an affine highest weight 



crystal (see Proposition 2.9) by analogy with the Kyoto path model, but now there are several 
highest weight components in the image, instead of just one. For each of these components we 
exhibit an explicit path from its highest weight (or ground state) to "type A elements" in the 
component, using only Demazure arrows (see Theorem |5.2[ ). This additional result suffices to 
establish the equality of energy and charge in type C, based on the corresponding result in type 
A. As a by-product, we obtain an explicit description of the components that appear in the 
nonperfect setting of single columns for type C. 

Our main result can now be stated as follows. 

Theorem 1.1. Let B = B"^^'^ ® ■ ■ ■ B^^'^ he a tensor product of KR crystals in type A^_-^ or 

type Cn with r^ > ttv-i > • • • > n > 1- Then for all b & B we have 

(1.3) D{b) = -charge(6), 

where D{b) is given in Definition \2A\ and charge(6) is defined in Section^ 

From a theoretical point of view, the above result is not surprising due to work of Ion [T2] , 
which relates Macdonald polynomials at t = and affine Demazure characters in simply-laced 
types. However, this result does not work in type C; in addition, it only gives an equality 
of polynomials (the generating functions for the statistics and the weights), not of individual 
terms. 

To compare our work with the previous papers on charge and the energy, let us first say 
that our results apply to arbitrary vertices in a tensor product of KN columns, not just to 
the highest weight elements (with respect to the nonzero arrows), that are used in the work 
involving Kostka-Foulkes polynomials. In type A, our approach via affine Demazure crystals 
comes closest to [19l HSJ 36]. However, in addition to the aspect mentioned above, it differs 
from the approach in these papers because we do not use the cyclage operation, which is based 
on the corresponding plactic relations; see [23] . These relations are the main cause of the 
complications in type C, in the work of Lecouvey |241 125| . 

The paper is organized as follows. In Section [2] we review the necessary crystal theory and 
define the energy function. In Section [3] we give the definition of charge both in types A and 
C. The proof of Theorem |1.1| for type A using the method of Demazure arrows is given in 
Section [4] In Section [5] we classify the various components under the Demazure arrows in type 
C in the extension of the Kyoto path model, and show that each ground state is connected 
to a type A filling. These results are used in Section [6] to prove Theorem 1.1 for type C. We 



conclude in Section [7] with a discussion of various possible extensions of this work. 
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2. Crystals and energy function 
In this section we review and set up crystal theory and define the energy function. 



2.1. Crystal generalities. Crystal bases provide a combinatorial method to study represen- 
tations of quantum algebras Uq{Q). For a good review on crystal base theory see the book by 
Hong and Kang [9|. Here g is a Lie algebra or affine Kac-Moody Lie algebra with index set /, 
weight lattice P, and simple roots Oj with i G I. The set of dominant weights is denoted by 
P+. For affine Kac-Moody (resp. finite Lie) algebras we denoted the fundamental weights by 
Aj (resp. Ui) for i £ I. 

A Q-crystal is a nonempty set B together with maps Ci, fi : B ^ B L) {0} for i £ I and 
wt: B ^ P. For 6 G B, we set ei{b) = max{k \ ef (6) / 0}, ipi{b) = max{k \ ff{h) / 0}, 

e{h) = ^£i{h)K and (/?(&) = ^fi{h)Ki. 

The beauty about crystal theory is that it is well-behaved with respect to taking tensor products. 
Let Bi and B2 be two g-crystals. As a set Bi <^ B2 is the Cartesian product of the two sets. 
For 6 = 61 (8) 62 £ -Si <8> B2, the weight function is simply wt(6) = wt(bi) + wt(62). The crystal 
operators are given by 

/.(6i0 62) = |^^('^^^^^ ifs.(M>v..(M, 
I 61 (X" fi{b2) otherwise, 

and similarly for ej(6). This rule can also be expressed combinatorially by the signature rule. 

A highest weight crystal B{X) of highest weight A G P-f is a crystal with a unique element 
ux such that ei{ux) = for all i G / and wt{ux) = A. On finite-dimensional highest weight 
crystals B{\) there exists an involution S : B{X) — )■ -B(A), called the Lusztig involution, which 
is a crystal isomorphism such that 

S{fi) = ei* and S{ei) = fi*. 

Here i* is defined through the map Oi 1— t- ai* := —wo{ai) with wq the longest element in the 
Weyl group of g. Explicitly, we have i* = n — z for type An-i and i* = i for type Cn- Under S 
the highest weight element goes to the lowest weight element. 

In [8], Henriques and Kamnitzer defined a crystal commutor on the tensor product of two 
classically highest weight crystals in terms of Lusztig's involution as 

B(\) (^ B{pi) ^ B{^l) ® B{X) 

^ ' ' bl(^b2^S{S{b2)®S{hi)). 
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2.2. Kashi^vara Nakashima columns for type C. Kashiwara and Nakashima [T7] de- 
veloped general tableaux models for finite-dimensional highest weight crystals for all non- 
exceptional classical Lie algebras g. For type C„, the Kashiwara-Nakashima (KN) columns 
|17j of height k index the vertices of the fundamental representation V{LOk) of the symplectic 
algebra sp2n(C^)- These columns are filled with entries in[n]:={l<2<---<n<n<n — 1< 
•••<T}. 

Definition 2.1. A column-strict filling b = 6(1) . . . b{k) with entries in \n] is a KN column if 
there is no pair {z,z) of letters in b such that: 

z = b{p) , z = b{q) , q — p < k — z . 

We will need a different definition of KN columns, which was proved to be equivalent to the 
one above in 



Definition 2.2. Let b be a column and I = {zi > ■ ■ ■ > z^} the set of unbarred letters z such 
that the pair {z,z) occurs in b. The column b can be split when there exists a set of r unbarred 
letters J = {ti > ■ ■ ■ > tr} C [n] such that: 

• ti is the greatest letter in [n] satisfying: ti < zi, ti ^ b, and ti b, 

• for i = 2, . . . ,r, the letter ti is the greatest one in [n] satisfying ti < min(tj_i, Zi), ti b, 
and ti b. 

In this case we write: 

• b^ for the column obtained by changing zi into ti in b for each letter Zi £ I, and by 
reordering if necessary, 

• b^ for the column obtained by changing Zi into ti in b for each letter Zi G /, and by 
reordering if necessary. 

The pair b b will be called a split column. 

Example 2.3. The following is a KN column of height 5 in type Cn for n > 5, together with 
the corresponding split column: 



L,R 



b'^b 



1 


4 


2 


5 


5 


3 


4 


2 


3 


1 



We used the fact that {zi > Z2} = {5 > 4}, so {ti > ^2} = {2 > 1}. 
We will consider Definition 12.21 as the definition of KN columns. 



2.3. Kirillov— Reshetikhin crystals. For the definition of the crystal energy function, we 
need to endow the KN columns with an affine crystal structure. These finite-dimensional 
affine crystals are called Kirillov-Reshetikhin (KR) crystals. Combinatorial models for all non- 
exceptional types were provided in [5]. 

Here we only describe the KR crystals B"^'^ for types ^„_i and Cn , where r G {1, 2, . . . , n— 1} 
and r G {1,2,..., n}, respectively. As a classical type An-i (resp. C„) crystal, the KR crystal 
is isomorphic to 
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The crystal operator /q is given as follows. Let b G B^'^, represented by a one-column 
KN tableau. In type A„_i, if b contains the letter n and no 1, fo{b) is the obtained from 
b by removing n and adding 1 to the column, leaving all letters in strictly increasing order. 
Otherwise /o(b) = 0. In type Cn, if b contains the letter 1, then fo{b) is obtained from b by 
removing the 1 and adding the letter 1, arranging all letters again in strictly increasing order. 
Otherwise fo{b) = 0. Note that if b contains 1, then it cannot contain 1 by the KN condition 
of Definition O 



Similarly, in type j4„_i, eo(6) changes a 1 into n if 1 is in 6 and n is not, and otherwise 
eo(6) = 0. In type Cn, eo(b) is obtained from b by changing a 1 into a 1 if it exists, and 
Co (6) = otherwise. 

2.4. The D function. Now we come to the definition of the energy function D on tensor 
products of KR crystals B^'^ of type ^^_i or Cn • It is defined by summing up combinatorially 
defined "local" energy contributions using the combinatorial i?-matrix. 

Let Bi , i?2 be two affine crystals with generators vi and V2 , respectively, such that Bi B2 
is connected and vi (X" V2 lies in a one-dimensional weight space. By [28', Proposition 3.8], this 
holds for any two KR crystals. The generator Vr,s for the KR crystal B^'^ is the unique element 
of classical weight scor- 

The combinatorial R-matrix [14, Section 4] is the unique crystal isomorphism 

(T : B2^Bi^ Bi® B2. 
By weight considerations, this must satisfy a{v2 ® vi) = vi (^ V2- 

As in |[14j and [39l Theorem 2.4], there is a function H = Hb2,Bi : -B2 SD -Bi — )• Z, unique up 
to a global additive constant, such that, for all 62 S B2 and 61 € i?i, 

— 1 if i = and LL, 

(2.2) H{ei{b2®bi)) = H{b2®bi) + {l if i = and RR, 

otherwise. 

Here LL (resp. RR) indicates that cq acts on the left (resp. right) tensor factor in both 
62 ® bi and a{b2 ® 61). When Bi and B2 are KR crystals, we normalize Hb2,Bi by requiring 
Hb2,Bx{v2 ® vi) = 0, where vi and V2 are the generators defined above. The map H is called 
the local energy function. 

Definition 2.4. For B = B''^'^ (g) ■ ■ ■ (g) B''^'^ of type A^^l-^ or ci^\ set 

Hji := ffjCrj+i(Tj+2 • • • crj_i and Hj^ := Hj^iaj^2<^j-3 • • • 0"j, 

where aj and Hj act on the j-th and {j + l)-st tensor factors. We define a right and left energy 
function D^, Dj^ : i? — )• Z as 

(2.3) D^:= Yl Hf^^ «^^ Db--= E ^f,^- 

N>j>i>l N>j>i>l 

We set Db := D^ and, when there is no confusion, we shorten Db to simply D; this is referred 
to as the energy function. 

Remark 2.5. Note that the energies D^ and Dj^ can be defined for general tensor products of 
KR crystals. When the KR crystals decompose into several components as classical crystals 
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(unlike in the cases for type ^„_i and Cn we consider), there is an extra term in the analogue 
of (2.3); see 



There is a precise relationship between D and D using the Lusztig involution. To state 
it, let us introduce the following map 

r : Bn ■ ■ ■ Bi ^ Bi ^ ■ ■ ■ Bn 

^AT (8) • • • 6i H> S{bi) • • • S{bN)- 

For types ^^_i and Cn and Bi = i?**"^, the KR crystal Bi is connected as a classical crystal 
and under S the classically highest weight element u-^^ '^^ maps to the lowest weight element 
^lowest j|. ^g j^q|. ]^g^j.(^ ^Q show from the explicit description of S, cq and /o in this case, that 
the following diagrams commute: 



(2.4) 



where B = B""^^^ 



s 

y 



eo 



5 

Y 



and 




(»5'^i'i andB = B''i'i 



^f-iv,!. 



This shows in particular that the crystal commutor (2.1) is lifted to an affine crystal isomor 



phism in these cases and hence must coincide with the combinatorial i?-matrix a. 
Proposition 2.6. Let B = B''^'^ (^ • • • (g) S^'i'^ of type A^l^ or C^^ and b£ B. Then 
(2.5) D§ib) = D^sirib)). 

Proof. Note that if a{bi (^ 62) = b'^ b[, then (7(5(62) S{bi)) = S{b[) 8(^2) since a and 



the commutor (2.1) agree on two tensor factors under the conditions of the proposition by the 



above arguments. Hence the terms in the sum of D^ are in one-to-one correspondence with 
terms in D^. Therefore it suffices to show that the local energy satisfies 



(2.6) 

for b2<E)bi £ B'^'^ 



Hib2 0bi) = HiSibi)^Sib2)) 



component as u 



highest 



B'-^'\ Since Siu^^^''') Si^'^^"''-^ - -1°-^* 



uT 



(g u. 



highest 



(g) U2°"'^^* lies in the same 



we have 



Hiu^'^^"'^ ^highest^ 



highest 



) = /7(5«s"'=^^) Siu^'^'"'''')) = 0. 



highest \ 



In the recursion (2.2), if for example cq acts LL on 62 (g bi, then /o acts RR on S{bi) ® £'(62) = 
S{bi) (8) 5(62) by (2.4). Hence H changes by the same amount in the left and right hand side 
of (2.6). The other cases can be checked analogously, which proves (2.6). D 



2.5. D energy as afRne grading. As suggested in \38\ Section 2.5] and proven in [32], the 
energy D is the same as the affine degree grading in the associated highest weight affine 
crystals up to an overall shift. We will explain this now since it plays a crucial role in the proof 
of the equality between charge and energy. 

We begin with the definition of Demazure arrows. For this we need constants Cr for r G -/^\{0} 
as for example defined in [Hj. In the cases of concern to us here, we have c^ = 1 for all r in type 
2 for 1 < r < n and Cr, 



Al^\ and c, 



1 in type Cn ■ 
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Definition 2.7. Let B = B"^^'^^ (^ • • • (8> B^"^'^'^ be a tensor product of KR crystals and fix an 
integer i such that i > [s^/cfc] for all 1 < k < N . We call such a tensor product a composite 
KR crystal of level bounded by £. 

An arrow fi is called an i-Demazure arrow on b £ B if ipi{b) > and either i £ I \ {0} or 
i = and So{b) > L 

In the setting of this paper, we are only concerned with tensor products of types ^„_i and 

Cn of the form B = B^^'^ (X" • • • (8) B^^'^. In this case one can pick 1=1 and a Demazure arrow 
for i? is a 1-Demazure arrow. 

Lemma 2.8. Let B = B''^'^ • • • (» B''^'^ of type A^^l-^ or cj^'^ and beB. Then 

(1) eo(6) > 1 implies D^{fo{b)) = D^{b) + 1; 

(2) ipo{b) > 1 implies D^{eo{b)) = D^{b) + 1. 

Proof. Part M follows directly from [42] Lemma 7.3]. For part Q, recall that by ( |2.4[ ) and 
Proposition |2.6| we have eo(T(6)) = t(/o(6)) and Z)^(r(6)) = D^{b). Also, setting b = r(6), we 
have ¥^0(6) > 1 if eo(^) ^ 1; thus, by using part ([l]), we deduce 

D\eo{b)) = D^ieoirib))) = D^(t(/o(6))) = D^(/o(6)) = D^{b) + 1 = D^ib) + 1 , 

which proves the claim. D 

The proof of the following essentially appeared in [14, Proof of Theorem 4.4.1] and was 
spelled out in this precise form in [42| Proposition 8.1]. Here 

P+ = {A G P+ I lev(A) = i}, 

where lev(A) := A(c) is the level of A and c is the central element c = J2iei '^t^t- 

Proposition 2.9. For B a composite KR crystal of level bounded by i, 

(2.7) B^B{iAo) = ^B{A'), 

A' 
where the sum is over a finite collection of (not necessarily distinct) A' G P^^ . 

In Section [5| we discuss in more detail the sum on the right hand side of (2.7) for type Cn 
and 1 = 1. 

Definition 2.10. For each b £ B, let u^ " be the unique element of B such that u^ ° '8' ii£Ao ^^ 
the highest weight in the component from Proposition 2.9 containing b(S>ui\g. 

Define the function deg on a direct sum of highest weight crystals to be the basic grading on 
each component, with all highest weight elements placed in degree 0. 

Corollary 2.11. Choose an isomorphism ra : B ® B{IAq) = @^i B{A'). Then for all b £ B, 
we have D{b) - D{ul^") = deg(m(6 (g) u^aJ). 

Corollary 2.12. The minimal number of eo in a string of et taking b to u^ " is D{b) — D[u^^ °). 

Remark 2.13. One special case of interest is when i? is a tensor product of perfect KR crystals 
of level (., and A also has level exactly i. Then the right side of (2.7) is a single highest weight 
crystal and the isomorphism in Proposition 2.9 is used i n the Ky oto p ath model [Hj. Hence 
u, " does not in fact depend on 6, simplifying Corollaries 



2.11 



and 



2.12 
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3. The charge construction 



3.1. The classical charge. Let us start by recalling the construction of the classical charge 
of a word due to Lascoux and Schiitzenberger [20]. Assume that u; is a word with letters in the 
alphabet [n] := {1, . . . ,n} which has partition content, i.e., the number of j's is greater than or 
equal to the number of j + I's, for each j = 1, . . . ,n — 1. The statistic charge(ii;) is calculated 
as a sum based on the following algorithm. Scan the word starting from its right end, and 
select the numbers 1,2,... in this order, up to the largest possible k. We always pick the first 
available entry j + 1 to the left of the previous entry j. Whenever there is no such entry, we 
pick the rightmost entry j + 1, so we start scanning the word from its right end once again; in 
this case, we also add k — j to the sum that computes charge(i(;). At the end of this process, we 
remove the selected numbers and repeat the whole procedure until the word becomes empty. 

Example 3.1. Consider the word w = 1132214323, where the first group of selected numbers 
is shown in bold. The corresponding contribution to the charge is 1. After removing the 
bold numbers and another round of selections (again shown in bold), we have 112323, so the 
contribution to the charge is 2. We are left with the word 123, whose contribution to the charge 
is 2 + 1 = 3. So charge(u') = 1 + 2 + (2 + 1) = 6. 

We now reinterpret the classical charge as a statistic on a tensor product of type Al^_^ KR 
crystals. Such a crystal indexed by a column of height k is traditionally denoted B^'^, and its 
vertices are indexed by increasing fillings of the mentioned column with integers in [n] . Given 
a partition fi (i.e., a dominant weight in the root system), let 



Ml 



(3.1) 



B,. 



B^^ 



i=l 



where /i' is the conjugate partition to fi. This is simply the set of column-strict fillings of the 
Young diagram /i with integers in [n]. Note that unlike in Section^ the tensor factors in (3.1) 
are ordered in weakly decreasing order. 



Fix a filling h in Bn written as a concatenation of columns 6i 



^/^i 



We attach to it a filling 



c := circ-ord(6) = ci 



-^1 



according to the following algorithm, which is based on the circular 



order -<i on [n] starting at i, namely i ^j i + 1 ^j • • • -<» n ^i 1 -<i • • • ^i i — 1. 

Algorithm 3.2. 

let c\ := h\; 

for j from 2 to fii do 

for i from 1 to fi', do 

let Cj{i) := min {hj \ {cj(l), , 

end do; 
end do; 
return c := ci . . . c^^. 



..,Cj(i-l)}, <cj_^{i)) 



Example 3.3. Algorithm |3.2| constructs the filling c from the filling h below. The bold entries 
in c are only relevant in Example |3.5| below. 



(3.2) 



3 


2 


1 


2 


5 


3 


2 




6 


4 


4 





and 



3 


3 


4 


2 


5 


2 


2 




6 


4 


1 
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We introduce some terminology in order to reinterpret the classical charge in terms of a 
statistic on B^. Given the considered filling b in S^, we define its charge word as the biword 
cw(6) containing a biletter ( .) for each entry k in the column bj of b. We order the biletters in 
the decreasing order of the fc's, and for equal fc's, in the decreasing order of j's. The obtained 



word formed by the lower letters j will be denoted by cw2(6). We refer to Example 3.5 for an 
illustration of the charge word. On the other hand, given the filling c = ci . . . c^^ constructed 



by Algorithm 3.2, we say that the cell 7 in column Cj and row i is a descent if Cj{i) > Cj+i(i), 
assuming that Cj+i{i) is defined. Let Des(c) denote the set of descents in c. As usual, we define 
the arm length arm(7) of a cell 7 as the number of cells to its right. 

It is not hard to see that Algorithm |3 . 2| for constructing c from b translates precisely into the 
selection algorithm which computes charge(cw2(6)). More precisely, consider the ith sequence 
1,2,... extracted from cw2(6) (which turns out to have length //j), and the letter j in this 
sequence; then the top letter paired with the mentioned letter j in cw{b) is precisely the entry 
Cj{i) in row i and column j of the filling c. In particular, the steps to the right in the ith 
iteration of the charge computation correspond precisely to the descents in the ith row of c, 
while the corresponding charge contributions and arm lengths coincide. We conclude that 

(3.3) y^ arm(7) = charge(cw2(6)) . 

7SDes(c) 

For simplicity, we set charge(6) := charge(cw2(6)). 

Remark 3.4. In [3^ we showed that the charge statistic on Bf^ can be derived from the Ram- 
Yip formula |40j for the corresponding Macdonald polynomial at t = 0. In fact, we showed that 



Algorithm 3.2 is closely related to the corresponding quantum Bruhat graph (see, e.g., p]). So 
we can conclude that this graph explains the charge construction itself. The mentioned idea 
was extended to type C, and it led to the definition of a type C charge, that we describe in 
Section 13.21 



Example 3.5. Note that cw2{b) for b in Example 3.3 is precisely the word w in Example 3.1 
In fact, the full biword cw(6) is shown below, using the order on the biletters specified above 
The index attached to a lower letter is the number of the iteration in which the given letter is 
selected in the process of computing charge(6): 

./654433222 1 
cw[b) -^^^ ^^ 3^ 23 2i li 4i 32 22 33 

One can note the parallel between the mentioned selection process and the construction of c 



from b in Example 3.3 The entries in the cells of Des(c) are shown in bold in (3.2). 



3.2. The type C charge. In this section we recall from [30| the construction of the type C 
charge. We start by fixing a dominant weight // in the root system of type C„. Let 

Ml 
(3.4) S^:=(g)5^^S 

i=l 

where B^'^ is the type Cn KR crystal indexed by a column of height k. Note that B^ is the set 
of fillings b = bi . . . b^-^ of the shape /i with integers in [n] whose columns bj are KN columns; 
indeed, the KN columns of height k label the vertices of B -^ . As mentioned above, it will be 
more useful to represent bj in the split form b^b^] in this case, b becomes a filling 6{'6f . . . b^^b^^ 
of the shape 2//. 
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Now fix a filling b in B^, represented with split columns, which are labeled from left to right 
by 1, 1', 2, 2', . • •• We can apply a slight modification of Algorithm 3.2 to b and obtain a filling 



c = CiC^ . . . cjl^cj}_^ = circ-ord(6) of 2fj,; namely, we start by setting c^ := bf, and then consider 
the (doubled) columns of b from left to right. We use the circular order on [n] starting at various 
values i, which we still denote by -<i. 

Example 3.6. Consider the following tensor product of KN columns: 



5 




[3" 




1 


3 

2 

1 


«) 


4 





3 




L3_ 




\l\ 



This is represented with split columns as the following filling b of the shape 2fi = (6, 6, 6, 2): 



1 r 2 2' 3 3' 



5 


5 


2 


3 


1 


1 


3 


3 


4 


4 


2 


3 


2 


2 


3 


2 


3 


2 


1 


1 











where the top row consists of the column labels. The corresponding filling c is 



5 


5 


4 


4 


3 


2 


3 


3 


3 


2 


1 


1 


2 


2 


2 


3 


2 


3 


1 


1 











Define the charge word cw(6) of b by analogy with type A, as the biword containing a biletter 
( .) for each entry k in column j of 6; here j and k belong to the alphabets {1 < 1' < 2 < 2' < . . .} 
and [n], respectively. We order the biletters as in the type A case (in the decreasing order of 
the fe's, and for equal /c's, in the decreasing order of j's), and define cw2{b) in the same way (as 
the word formed by the lower letters j). 

The modification of Algorithm |3.2| for constructing c from b can be rephrased in terms of 
cw(6), as explained below; we will refer to this rephrasing as the charge algorithm. We start 
by scanning cw2(6) from right to left and by selecting the entries 1, 1', 2,2' , . . . , fii, (/^i)' in this 
order, according to the following rule: always pick the first available entry to the left, but if the 
desired entry is not available then scan the word from its right end once again. As in type A, 
we can see that the sequence of top letters paired with 1, 1', 2,2' , . . . , fii, (/xi)' is the first row 
of the filling c (read from right to left). We then remove the selected entries from cw2(6) and 
repeat the above procedure, which will now give the other rows of c, from top to bottom. It 
was shown in [30] that we always go left from j to j' , but we can go right from j' to j + 1. 



Example 3.7. This is a continuation of Example 3.6 The charge word cw(6), with the order 
on the biletters indicated above, is 



1 

1'. 



1222233334455332211 
1-4 I4 3']^ 22 I3 I3 3i 22 I2 I2 2'i 2\ V^ \\ 33 23 33 23 32 32 

The index attached to a lower letter is the number of the iteration in which the given letter is 
selected by the charge algorithm. 



Descents are defined as usual, cf. Section [37Lj It is easy to see that the descents in c correspond 
to the steps to the right in the charge algorithm applied to cw2(&). By an observation made 
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above, we only have descents of the form c^{i) > Cj^i{i). We are led to the following definition 
of the type C charge. 

Definition 3.8. Consider a word w with letters in the alphabet 1, l', 2, 2', . . ., containing as 
many letters j as f , and at least as many letters j as j + 1. Apply the charge algorithm to 
w, and assume that a selected entry j' is always to the left of the previously selected j. Let 
charge(t(;) he the sum of k — j for each selected entry j + 1 to the right of the previously selected 
f , where the selected entries in the given iteration are 1, 1', . • • ) ^; k' . 

The above discussion leads to the following result: 

(3.5) - y^ arm(7) = charge(cw2(6)) . 

7£Des(c) 

For simplicity, we again set charge(5) := charge(cw2(6)). 



Example 3.9. This is still a continuation of Example 3.6 The entries in the descents of c 
are shown above in bold. Correspondingly, the charge algorithm applied to cw2(6) makes one 
step to the right in the second iteration (from 2' to 3), and two steps to the right in the third 
iteration (from 1' to 2 and from 2' to 3). Thus, charge(6) = 1 + (2 + 1) = 4. 

4. Energy and charge in type A 

In this section we rederive the result of Nakayashiki-Yamada [37] showing the equality of 
the energy function and charge in type An-i- We do this in a more conceptual way, by using 
the method of Demazure arrows (the proof in |37j is based on subtle combinatorics of Young 
tableaux). Furthermore, we work with all the crystal vertices in a tensor product of columns, 
not just the highest weight vertices considered in [37]. Note that the setup in the mentioned 
paper is that of the right energy and a tensor product of columns of increasing heights which, 
by Proposition |2.6[ is equivalent to the setup in this paper. 



We start by studying the behavior of the type A charge with respect to the crystal operators. 
Proposition 4.1. The type An~\ charge is preserved by the crystal operators /i, . . . , fn-i- 



Proof. Let 6 be a tensor product of columns in some S^ (see Section 3.1 and the terminology 

therein, which we use freely). It is known that the word cw2(/i(6)) is in the same plactic 

equivalence class (see, e.g., [21j ) as cw2(6). More precisely, the former is obtained from the 

f i + 1 
latter by considering its subword formed by the letters x corresponding to biletters 



X 

or I ] in cw(6), by viewing this subword as (the column word of) a skew tableau with two 

columns, and by using jeu de taquin to move a letter from the right column to the left one. 
This is explained in detail in |22i Section 2], based on the notion of "double crystal graphs". 
Then we use the well-known fact that the classical charge is preserved by the plactic relations 
(see, e.g., [211 Lemma 6.6.6 (ii)]). □ 

Proposition 4.2. Let B = B^'^'^ ® • • • (gi B^^'^ be of type ^„_i with r^ > r^-i > • • • > ?'i and 
b & B. If ipQ{b) > 1 and £o{b) > 1, then the type An-i charge satisfies 

charge(eo(6)) = charge(5) — 1 . 
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Proof. Let b = hi . . . h^^ , and assume that cq changes the entry 1 in bj to n. The condition 
Vo{b) > 1 imphes j > 1. Let c := circ-ord(6) with c = ci...c^^, where Cj{i) = 1, and 



d := circ-ord(eo(6)) with d = di . . . d^-^ (recall that the map circ-ord is defined by Algorithm 3.2). 

The main case. We start by assuming that Cj-i{i) > 1 and that the column Cj+i, if it exists, 
does not contain n in a row k > i. We claim that, in this case, di{k) = ci{k) for all k,l with 
the exception of dj{i) = n. Indeed, all we need to check are the following, in this order: (1) the 
entries dj{k) for A; = 1, . . . , i — 1 are as claimed, since the alternative, namely that one of them 
is n, would lead to a contradiction; (2) the value of dj{i) follows from the fact that Cj-i{i) / 1; 
(3) the value of dj+i{i), if this entry exists, follows from the above condition on the column 
Cj-|_i. By the above facts, we have the same descents in the fillings c and d, with the exception 
of the descent Cj-i{i) > Cj{i) = 1, which corresponds to dj{i) = n > dj+i{i), assuming that 
dj+i{i) exists. The difference between the arm lengths of the mentioned descents is 1, which 



concludes the proof by (3.3). The remaining part of the proof treats the exceptions to this case. 



Exception 1. Assume that Cj_i(i) = 1. Then the column bj-i must contain n (otherwise the 
entry 1 in bj would not be the leftmost unpaired 1). Assuming that Cj-i{k) = n, we must have 



fc > i, by Algorithm 3.2 and the fact that Cj{i) = 1. The condition (fo{b) > 1 implies j > 2, and 
we have Cj-2{i) = 1 (otherwise Cj_i(i) = 1 and Cj-i{k) = n for k > i are in contradiction with 



the way in which Algorithm 3.2 reorders the column fcj-i). This reasoning can be continued 



indefinitely, so this case is impossible. 

Exception 2. Assume that Cj-\{i) > 1, but Cj+i{k) = n for some k > i. Then bj^i must 
contain 1 as well, namely Cj+i(Z) = 1 (otherwise the entry 1 in bj would be paired with the 



entry n in 6j+i). We must have / < i, by Algorithm 3.2 and the fact that Cj{i) = 1. The above 
facts imply that I < k. 

Exception 2.1. Assume that I < i. The fact that Cj{l) > 1, Cj+i{l) = 1, and Cj+i{k) = n for 
k > I are in contradiction with the way in which Algorithm |3.2| reorders the column 6j+i. So 
this case is impossible. 

Exception 2.2. The only possibility left is that Cj+i{i) = 1 and Cj+i{k) = n for k = ki > i. 
Let us assume for the moment that the second condition in Exception 2 holds for the column 
Cj_|_2, namely 0^+2(^2) = n for some k2 > i- By the same reasoning as above, we deduce that 



Cj_|_2(i) = 1, so ^2 > i. In fact, we also have fci > k2, by Algorithm 3.2 By continuing this 
reasoning, we obtain 



Cj+i{i) = ... = Cj+p{i) = 1 , Cj+i(A;i) = . . . = Cj^p{kp) = n , where ki> ...>k. 



•p ) 



on the other hand, we can assume that, if the column Cj+p+i exists, then it does not contain 
n in rows i,i + 1, . . .. This information about the filling c is represented in the figure below. 
The column Cj is the leftmost column with an entry 1 displayed, while + and * stand for 
entries different from 1 and n, respectively. The boxes shown in bold represent descents. By a 
reasoning similar to the main case above, we deduce that the only difference between the fillings 
c and d consists of the entries 1 and n in the figure changing to n and 1 in d, respectively. This 
leads to the marked descents moving to the boxes indicated by the arrows. It is now easy to 
see that the sum of the arm lengths of descents decreases by 1 when passing from c to d. 
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1 1 



+ n 



1 




n 


* 



+ n 



1 
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* 



+ n 



1 


* 






n 


* 



D 



tions 



4.1 



4.2 



Proo f of Theo rem 1.1 in type ^4. The proof is immediate based on Corollary 2.11 and Proposi- 

D 



(1) 



using the fact that KR crystals of type A^_^ are perfect. 



5. Kyoto path model for nonperfect type C 



.(1) 



In this section, we make Proposition 2.9 more explicit in the case of i? = B'''^'^ 



^n 



and £ = 1 for type Cn , by providing a correspondence between highest weight elements (or 



ground states) in B (^ B{Aq) and elements in B{A') in the sum on the rig ht h and side of (2.7), 
which are of type A. This will help in the next section to prove Theorem |l.l| for type Cn ■ 

We call the highest weight elements in i? (8) -B(Ao) ground state paths. There is a recursive 
construction for them, which starts by listing all elements bi £ B^^'^ such that e(6i) = Aq. 
Suppose 6fc (8) • • • (8) 6i G B^'"^ (8 • • • (8 B'^^'^ are already constructed. Then bk+i G B^'=+^'^ can be 
any of the elements such that e{bk+i) = ^{bk)- The weight of the ground state is if{biy), which 
is some fundamental weight Ah- For perfect crystals there are unique elements b]\f, . . . ,bi with 



.(1) 



the described properties. However, in type Cn the crystals B^' are not perfect and the above 
construction gives a tree of ground state elements. 



Example 5.1. Take B = B 



1,1 



B 



2,1 



^2,1^^3,1 of type d^\ 



Then bi is the column 321 



and 62 the column 23. For 63 there are two choices, namely the columns 32 or 22. In the first 
case 64 is 3, and in the second case 64 can be 2 or 1. In summary the three ground states are 



3 ® 


2 
3 





3 

2 





1 
2 
3 



^Ao 



s 



« 


2 


CR) 


3 


(8 


1 

2 




2 




2 




3 



The weights are A2, A2, and Aq, respectively. 



(8nAo 



.(1) 



1 (8 


2 


^ 


3 


CR) 


1 




il_ 




2 




3 



(8U.A0 • 



Theorem 5.2. Let B = B^^'^ (8 • • • <8 B^^'^ of type Cn ■ From each ground state u (8 u\q G 
B (8 -B(Ao) there exists a sequence of Demazure arrows fi (see Definition 2.7), which ends at 
an element b (8 "UAq such that b does not contain any barred letter. 



Note that in Theorem 5.2 there is no assumption on the order of the r^. 
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In order to provide a proof of Theorem 5.2, we describe the exphcit sequence of fi satisfying 
the required conditions. For this we recursively define the following objects. 

- Let Ao ^ Ai C • • • C As (s is given below) be a sequence of shapes, where Aq is a single 
column of height h if the weight of the ground state path u <^ u/\_^ is A/j. The other 
shapes are all of the form A(A;, /12, /ii), a partition with k columns of height n followed 
by two columns of heights n > /i2 > /ii > 0, respectively. If Xj = A(A;, /12, /ii), then 
Aj+i = A(fc, /i2 + 1, /ii + 1), where we identify X{k, /i2 + 1, /ii + 1) = A(A; + 1, /ii + 1, 0) 
if /i2 + 1 = n and /ii + 1 < n, and \{k + 2,0,0) if h2 + 1 = hi + 1 = n. This adds one 
horizontal domino in the consecutive rightmost columns up to height n. 

- We recursively define 

Vj+i ■.= Fj{vj) 
with vq = u the ground state and Fj the sequence 



fc+l r2k+2 



] •~ 



Jo /l 



2k+2 f2k+l 
hi+1 



fZK+Z f 

Jhi Jl 



2k+l f2k 



Jh2 Ji 



e2k 



h2 + l 



Jn—lJn 



if Xj = X{k,h2,hi). We continue doing this as long as possible, in other words, until 
FJv^) is undefined. 



Example 5.3. Let u uaq be the second ground state from Example 5.1 Then the sequence 
of Ao C Ai C • • • C As is 



C 



c 



c 



with 



/0/1/2 



2 f3 f2 . 



Fi — /o /i /2 h 



2 f4f3 . 



F2 — /o /l /2 h 



Let Xj = (Aji > Aj2 >■■■), and consider the conjugate partition A'- = (A'-j^ > A'-g >•••)• We 
will now explain how to represent an element Vj by a collection of non-crossing lattice paths 
(which might touch each other); see Figure [Tl The paths have the following types of steps: down 
(southeast), up (northeast), and horizontal (east). More precisely, the paths in the collection 
Pj = {Pi) • • • J Pi} representing Vj correspond to the / = Xji columns of Xj, and they satisfy the 
conditions below. 



(1) The endpoints oi pi, . . . ,pi are aligned at height from right to left, and pi starts at 
height X'j^. 

(2) The paths p2, ■ ■ ■ ,Pi and the segment of pi before the end of p2 (in case p2 exists) 
consist entirely of down and horizontal steps. Two up or down steps never lie below 
one another, and neither do only horizontal steps. 

(3) The path pi starts after pj+2 ends, for i = 1, ...,/ — 2. 

The terms "before" and "after" refer here to the order of the corresponding x-coordinates, with 
equality allowed. Each down (resp. up) step starting at height i is labeled by a letter i (resp. 
i). We define word(Pj) as the word obtained by reading the labels on the paths in Pj from left 
to right (recall that two labels never lie below one another). 

Let us now explain the construction of the collections Pj . We start by defining Pq as consisting 
of a single path: the one having the same word as u = fo, ending at height 0, and having no 
horizontal steps; the word of u, denoted word(u), is defined as usual, by reading its columns 
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from left to right, bottom to top. We construct the collections Pj recursively, via the following 
transformation rule P i— )■ up(-P) on a collection of paths P = {pi, . . . ,pi) satisfying conditions 
(l)-(3) above. 

Rule 5.4. 

- Locate the leftmost up step inpi, and let y he its label. Replace it with a horizontal step 
at height y and a down step below it ending at height 0. 

- Shift up by 1 the segment ofpi to the left of the position where the above change occurred, 
to connect it to the tail of pi. Shift up by 1 all the other paths. 

- Replace any down step above height n by a down step below it, ending at height 0. 

- Consider the down steps ending at height from right to left, excluding the rightmost 
one. Match them with the shifts of p2, ■ . . ,Pi, in this order, and connect the matched 
pairs by horizontal lines of height 1 . 

Note that, in the last step of the rule, the last one or two down steps might have no match, 
so they start new paths. For simplicity, any horizontal steps at the beginning of a path are 
ignored. It is easy to see that, since P satisfies conditions (l)-(3) above, the rule can be applied, 
and up(P) satisfies the same conditions. Thus, we can recursively define -Pj+i := up(Pj) as 
long as there are up steps in Pj. 

We claim that applying the rule to Pj corresponds to applying Fj to Vj. To make this precise, 
we introduce some notation. Let r be the collection of columns oi B = B"^^'^ • • • B^^'^, i.e., 
the tuple (rTv, . . . , ri). Given a word w of length ri + • • • + r^v, we define a rearrangement of 
it ord,.(u') by slicing w into segments of length r^r, . . . , ri in this order, and by reordering each 
segment decreasingly. We will show below that, for all j = 0, . . . , s, we have 

(5.1) ordr(word(Pj)) = word(t'j) . 

We will see that the reordering of the segments of word(Pj) done by ord,-, if any, is very 
simple: a segment is a concatenation of two decreasing segments which can be swapped to give 
a decreasing sequence. Based on the above discussion, one can see that Xj is the weight of Vj, 
and that this element is of highest weight in its classical (non-affine) component. 

Example 5.5. The collections of paths Pj corresponding to Example 5.3 is shown in Figure [T] 



Proof of Theorem 5.2 . It essentially suffices to show that Fj (vj) is defined if and only if up(Pj) is 



defined and that ( 5.1 ) holds, for all j (the only extra fact to check is that all arrows corresponding 



to Fj are Demazure arrows, see below). Indeed, Rule 5.4 can be applied as long as there are 
up steps in Pj, so the final element Vg will have only positive entries. We will first show that 
Rule 5.4 precisely describes the action of Fj on the word of Pj, viewed as an element of the 



tensor product (^B^'^)'^^^^~^ i-''iv). more precisely, we have 

(5.2) Fj{woid{Pj)) = word(Pj+i) . 

We will then show that the above action of Fj is completely similar to that on Vj, which will 



prove (5.1). 



Let us describe the action of Fj on word(Pj). First note that the connected components of the 
paths in Pj correspond to bracketed units in the tensor product (i?i'i)®(^'i"* ^''''^l according to 
the signature rule for the application of the Kashiwara operator on tensor products of crystals. 
This means that any given operator /j can only operate on the starting points of the connected 
paths in Pj. The application of f^ lifts all starting points of value n to n. Note that by the 
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Figure 1. Paths for Example 5.5 



weight \j = X{k, /i2, hi) of Vj there are precisely k of them. Next /^^i lifts all A; just created 
n to n — 1 and the k letters n — 1 to the left of the leftmost up-step to n that were previously 
bracketed with the just lifted n etc.. Potentially there are two more connected path components 
with starting points at hi and /12, hence the exponents of fi with i < hi and i < h2 are increased 
accordingly. In summary, the starting points n of path components are eventually lifted to 1 
and then turned into 1 by /o; the other positive letters to the left of the leftmost up step are all 
raised by one. In the rightmost path, this lifting process eventually reaches the leftmost up step 
y which is lifted to y — 1 by fy-i- The steps to the right of the leftmost up step y are either 
y + 1 or y. Both would need an fy to be lifted. Since the fi in Fj are applied in decreasing 



order of the indices i, this cannot happen and hence the leftmost up step is lifted to 1 and then 



turned into 1 by /q. Altogether, the changes are precisely as described in Rule 5.4 
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Let US now compare the action of Fj on word(Pj), i.e., on (i?l■l)®(''l^ t-'^iv)^ with that on v 



J' 



i.e., on i? = B^'^'^ (S" • • • (8> -B''^'^. Note that by the action of /o on columns rather than tensor 
products of single boxes, the word word(Pj) of a path is a concatenation of cyclically shifted 



columns of Vj. 



We claim that Fj acts on precisely the same entries in the two cases. In the 



expression for Fj the /« with larger indices i act first and hence lift up the upper or barred 
portion of the paths. This is the case before or after reordering, /o changes a 1 into a 1 in both 
cases. Hence the action commutes with the reordering. 

Finally, the rightmost tensor factor in vq is the column c = ri • • • 1 G B"^^'^. Note that this 
column is never changed during the algorithm and satisfies eo(c) = 1. By the tensor product 
rules this implies that eo{vj) > 1 for all < j < s, so that all arrows are Demazure arrows. D 



There is a more direct way of constructing the type A elements b in Theorem 5.2 from the 
ground state u(^u/\_g. For i = 1,2, . . . ,N, place the letter i in column 1 for each unbarred letter 
in Ui and in column 2 for each barred letter, where u = un (^ ■ ■ ■ (^ ui. Note that, due to the 
fact that u is a ground state, the difference in height between the first and second column is 
at most n. Now cut the columns at heights n, 2n, 3n, . . . and put all pieces next to each other 
aligned at height 0. The letters in row i record the tensor factors of b in Theorem 5.2 which 
contain the letter i. Comparing this with the algorithm in the proof of Theorem 5.2 
hard to see that it gives the correct answer. 



it is not 



Example 5.6. Continuing Examples 5.3 and 5.5, the chosen ground state yields the columns 



4 
3 




1 


3 


1 


2 


1 


2 



1 


3 




1 


2 


4 


1 


2 


3 



and, after the cut, 



This tells us that the letter 1 appears in bi, b2, 63, the letter 2 appears in 61, ^2, ^4, and the 
letter 3 appear in 61 and 63. This agrees with Figure [T| 



6. Energy and charge in type C 

The purpose of this section is to provide the proof of Theorem 1.1 for type C. 

6.1. The energy for type A fillings. We start with a brief discussion of the combinatorial 
i?-matrix Bi®B2 -^ B2® Bi for a tensor product of two type C columns. By the type C Fieri 
rule [l9], the decomposition of Bi ® B2 into classical (non-affine) components is multiplicity 
free. On the other hand, it was proved in [23] that the type C jeu de taquin due to Sheats 
|44j is compatible with the classical crystal operators. We conclude that the mentioned jeu 
de taquin, when applied to our two-column situation, realizes the corresponding combinatorial 
i?-matrix. Note that, in all situations considered below, the type C jeu de taquin is governed 
by the simpler rules of the classical one, in type A (see, e.g., ^). 

Lemma 6.1. Let b = bi®b2 be a tensor product of two columns with all entries in [n], where 
the height 0/62 is at most that ofbi. Then the type ^^_i and Cn (local) energies ofb coincide. 

Proof. Consider the crystal B^'^ (8) i? 2.1 of type Cn , where n > /ii > /12 > 1. By the type 
C Fieri rule, the classical components of this crystal containing fillings with all entries in [n] 



have highest weights fii = {hi + i,h2 



for 



0, . . . ,min(/i; 



2,n 



/ii). We will calculate the 



CRYSTAL ENERGY VIA CHARGE 



19 



(local) energy on these components based on its definition (2.2), by setting the energy to on 
the component of highest weight /io (the usual normalization). 

For i = 1, . . . ,min(/i2,n — hi), consider bu b2i in B^^'^ B^"^'^ , where 



n — hi + 1 



bu 



n 



b2i 



1 


n-hi-i + 2 




n — hi 


n - /i2 + 1 




n — i 



The image bgj ® b[^ of bu (E) 62* under the combinatorial i?-matrix, constructed with jeu de 
taquin, is given by 



n - /i2 + 1 



b'. 



2i 



n 



b'u 



1 


n - 


-hi-i + 2 




n — i 



Clearly, cq acts RR on bu ® b2i and 62^ (E> b[-, by changing 1 to 1. By type C insertion [23] . 
which in this case is essentially just type A insertion (see, e.g., [7]), it is easy to see that bu®b2i 
and eo(6ii ® 62*) lie in the components of highest weights ^i and /^i-i, respectively. Therefore, 
the energy on the component of highest weight Hi is —i. This coincides with the type A energy 
(calculated via a similar procedure or via the type A charge). D 

Proposition 6.2. Ifb is a tensor product of columns with all entries in [n], then the type An-i 
and Cn energies of b coincide. Furthermore, if the columns have weakly decreasing heights, they 
equal the (type An-i or Cn) charge of b. 

Proof. When all entries are in [n], the jeu de taquin algorithms in types vln-i and C„ (realizing 
the corresponding combinatorial i?-matrices) work identically. Since the corresponding local 



energies coincide by Lemma 6.1 , the global energies coincide as well. But the type An-i energy 
of b is computed by the type An-i charge when the heights of the columns are weakly decreasing 
by Theorem 1.1 which was proven in Section |4j which clearly coincides with the corresponding 
type Cn charge. D 



6.2. The conclusion of the proof. We start by studying the behavior of the type C charge 
with respect to the crystal operators and state a result in [3T] . 

Proposition 6.3. [31j The type Cn charge is preserved by the crystal operators /i, . . . , /„. 

Proposition 6.4. Let B = B^'^'^ (8) • • • (8) B^^'^ be of type Cn with rjv > t'at-i > ■ ■ ■ > ri and 
b & B. If ipQ{b) > 1 and eoib) > 1, then the type Cn charge satisfies 

charge(eo(6)) = charge(5) — 1 . 

Proof. Let b = bi . . . 6^^, and note that none of these columns contain both 1 and 1. Assume 
that eo changes the entry 1 in bj to 1. The condition ipo{b) > 1 implies j > 1. Note that the 
column bj-i, and thus b^_^, cannot contain 1 (otherwise, the entry 1 in bj would not be the 
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leftmost unpaired 1). Similarly, the column 6j+i, and thus &£^i, if they exist, cannot contain 1 
(otherwise, the entry 1 in bj would be paired with the mentioned 1). 

Let c := circ-ord(6) with c = CiC^ . . . cj;^-i^cj}_^ , and d := circ-ord(eo(6)) with d = dfd^ . . . dj;^_^dj^_^ 



(recall that the map circ-ord is defined by a slight modification of Algorithm 3.2 ). Ass ume that 
Cj{i) = 1, which implies c^(i) = 1 (as we have no descents in the pair c'j'Cj^; see Section 3.2). We 
claim that df"{k) = c('{k) and df{k) = cf-{k) for all fc, I with the exception of d^{i) = df{i) = 1. 
Indeed, all we need to check are the following, in this order: (1) the entries d^{k) and d^{k) 
for k = 1, ... ,i — 1 are as claimed, since the alternative, namely that one of them is 1, would 
lead to a contradiction; (2) the value of d^{i) follows from the fact that b^_^ contains no 1; (3) 
the value of d^j^^{i), if this entry exists, follows from the fact that 6:L_i contains no 1. 

By the above facts, we have the same descents in the fillings c and d, with the exception 
of the descent c^_i{i) > c^{i) = 1, which corresponds to d^{i) = 1 > d'j'_^_i{i), assuming that 
df_^^i{i) exists. The difference between the arm lengths of the mentioned descents is 2, which 



concludes the proof by (3.5). D 



Proof of Theorem 1.1 in type C. The proof is immediate based on Lemma 2.8, Theorem 5.2 



and Propositions 6.2 6.3 6.4 □ 



7. Open problems 

In this section we discuss several directions of research stemming from the results in this 
paper. We intend to pursue some of them in the future. 

7.1. Columns of types B and D. We believe that, in types B^ and Dn, the statistic in the 
Ram- Yip formula for Macdonald polynomials at t = can also be translated into a charge 
statistic on the corresponding tensor product of KR crystals B'^'^ of types Bn and Dn (rep- 
resented with KN columns), thus extending the results in [30]. We conjecture that the type 
Bn and Dn charge also agrees with the corresponding energy function, which would extend 
the results in this paper. It should be possible to use a proof technique similar to the one 



in this paper, based on Lemma 2.8 from ^42j, to prove the conjecture; the proof would in- 
clude the generalization of the results in [3l] to types Bn and Dn- These claims are supported 
by |42t Corollary 9.5] expressing a type Dn Macdonald polynomial at t = in terms of the 
corresponding energy function. 

Prom one point of view, the case of types Bn and Z?„ is easier than the one of type C„, 
because the corresponding level 1 KR crystals are perfect. However, the construction of charge 
in types Bn and Dn displays additional complexity, due to some new aspects, that we now 
describe. We start by referring to type Bn, as type Dn has all the complexity of type Bn plus 
an additional one. Por the corresponding KN columns, indexing the vertices of the crystals 
B(ui) corresponding to the fundamental representations V{ijJi), we refer to [17j. 

The first new aspect in type Bn is related to the splitting of the KR crystal i?'^'^ upon 
removing the 0-arrows as the following direct sum of (classical) crystals: B{ijJk)®B{ujk^2)®- ■ ■ ■ 
This phenomenon manifests itself in the existence of descents between the left and the right 
columns of a split KN column, which was not the case in type Cn- We illustrate this based on 
the following example. 
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Example 7.1. 

column 



Let n = UJ4 in i?5. The KR crystal B^'^ contains a vertex indexed by the KN 



e B{uj2) C B*'' ~ B{ooi) e ^(wa) S(wo) • 



Instead of constructing the split column, we first construct from b the following "extended" KN 
column of height 4, in order to make the columns corresponding to the components of B"^'^ of 
the same height: 



this construction is based on a procedure in [HI Section 3.4]. Only at this point we construct 



the split column. We claim that the doubling procedure in Definition 2.2 can be extended, 
but we still match the entries in / = {3 > 1} with certain entries "preceding" them; the only 
difference is that this now requires us to go counterclockwise around the circle. So we obtain 
J = {2, 4}, which gives the following splitting of h by the usual rule: 



LtR 



h% 



2 


1 


4 


3 


3 


4 


1 


2 



Finally, we construct c 



c^c^ 



circ-ord(6'^6^) by the usual rule (see Section 3.2), where 



2 


3 


4 


2 


3 


1 


1 


4 



The energy of 6 in B^'^ is —1. This agrees with the charge oib b , computed via the sum in (3.5); 
indeed, there are two descents in c, whose arm lengths are 1, so charge(6^6'^) 



5(1 + 1) 



1. 



The second new aspect in type Bn is the fact that not always a descent in the filling obtained 
via the procedure "circ-ord" contributes half its arm length to the charge. We illustrate this 
with another example. 

Example 7.2. Let fi = 2oj2 in type B^, so Bn = B^'^ B'^'^. Consider the pair of KN columns 



1 


1 


2 


2 



in B^. The corresponding fillings with split columns b = 6f'6{^62^2^ ^'^d c = circ-ord (6) 
cf'cf C2 c^ coincide, and they are 



1 


1 


1 


1 


2 


2 


2 


2 



The corresponding energy is —2, but in c we have only one descent of arm length 2, so we 



cannot use the sum in (3.5) to compute the charge. This problem occurs because, for all the 
entries i between 2 and 2 in circular order (namely i G {1, 1}), either i ox i are above 2 in c{^. 
This forces the descent 2 > 2 to be double counted. A more illuminating explanation can be 
given in the context of deriving the charge from the statistic in the Ram- Yip formula; see [30^ 
Section 9]. 
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The above example suggests that in type Bn we need to modify the definition of charge by 



the sum in (3.5) as follows. We use the same notation c = CiCi ■ ■ ■ c^-i^c^_^ as above. Let Des'(c) 
denote the descents of the form m = c^{i) > Cj'j^i{i) = m such that, for any A: = 1, . . . , m — 1, 
we have either k or k in cf'[l,i — 1]. We claim that the appropriate definition of the type B^ 
charge is given by the following formula: 

(7.1) - ^ arm(7) + ^ arm(7) . 

7eDes(c)\Des'(c) 7eDos'(c) 

We said above that in type Dn we have additional complexity still. We explain this using 
another example. 

Example 7.3. Let /x = 2uji in type D^. The filling with split columns c = b = 6f'6{^62^|^ = 



3 3 3 3 in B^ has no descents, but the corresponding energy is —1. The reason for this is 



that the values 3 and 3 are incomparable in type D^, so the pair 3 3 in c needs to contribute 1 
to the charge. In fact, the definition of charge needs to be adjusted even more, as the following 
more subtle example shows. 

Let fj, = 2uj2 in D4. Consider the following filling with split columns in i?^: 



bfbfb^b^ = I 



3^1 
4 3 



This filling has only ascents or equal entries next to each other in a row. However, the corre- 
sponding energy is —1, so the charge needs to be 1. The reason for which the charge is not 
0, meaning that b is not a split KN tableau, is that the two middle columns in b represent 
a forbidden configuration in type D4. (Recall from [17] that, in type D, in addition to the 
row monotonicity condition for the split tableau in types B and C, there are extra conditions 
for a sequence of KN columns to form a KN tableau, and these are given in terms of certain 
forbidden configurations for a pair formed by a right column and the next left column.) 

7.2. Rows of types B, D, and C. In type A, Nakayashiki and Yamada [37] showed that the 
Lascoux-Schiitzenberger charge statistic is related to the energy function on tensor products of 
both rows and columns. In particular, in [37' Proposition 3.23] they provide an explicit relation 
between the energy function on rows and columns by giving a bijection of each set to the set 
of semistandard Young tableaux which preserves the statistics. 

We expect a similar result to hold in the other classical types (with slight modifications). 
This claim is motivated by |28| Theorem 10.10], which relates the one-dimensional configuration 
sum in the large rank limit for columns to the one for rows in the dual type. More explicitly, 
denote by X^"^{q) the one-dimensional configuration sum, that is, the sum of q^^^> over all 
highest weight elements b in the crystal B of weight A and type Yn, graded by the energy 
D{b). Then up to an overall power of q, X^'^{q) and Xy"^,{q~^) are equal for large rank n. 
Here, if B = B''^'"^ (g) ■ ■ ■ (g) B^^^''^ then B' = 'b"^^''^ (g) ■ ■ ■ ® B''^'"'^ and A' is the transpose of A 
(where we identify dominant weights with partitions). Note that X^'^{q) = X^^ 'j^{q) when n is 
large [281 [50]. 

I n pa rticular, if the charge for columns of types i?„ or D„ were known as outlined in Sec- 
tion 7.1 then, using i? = B^'^'^(g- ■ ■(gB'''^'^, the relation between X_)^^(g) and Xy"^,{q'^^) would 
provide a way to obtain the charge in type Cn also for single rows (as opposed to single columns 
as treated in this paper). Computer experiments using Sage [47J indicate that this duality not 
only holds on the level of configuration sums but also on the level of individual terms. We 
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expect that the highest weight crystal elements for both columns and rows are in bijection with 
certain oscillating tableaux, as introduced in [1], due to the fact that the recording tableaux for 
the Robinson-Schensted algorithm for types B^, Cn, D^ as introduced in [231 1221 are given by 
these tableaux. If this conjecture were to be true, then we would be able to obtain the charge 
for single rows from the charge for single columns via this bijection in the stable range (that is, 
large n). 

7.3. Arbitrary order of columns and arbitrary types. Consider a composition r = 

(rjv, . . . , ri), and the crystal B := B^'^'^ (d • • • B"^^'^ in types A and C. Since the charge 
is defined only if r is a partition, the computation via the charge of the energy function on 
some vertex b in B must proceed in general indirectly, by first appropriately commuting the 
columns of b via the combinatorial i2- matrix (which preserves the energy). However, it is de- 
sirable to compute the energy directly on i3. In addition, we would like to generalize to B of 
arbitrary type, i.e., to find a statistic on B expressing the energy, which can be computed using 
only the combinatorial data corresponding to a vertex, instead of various paths in B. 

Let fi{r) be the partition obtained by reordering the parts of r. The Ram- Yip formula for 
the Macdonald polynomial P^(r)(a;; g, 0) can be written in a way that is compatible with the 
composition r rather than the partition fi{r). Recall from [30t Proposition 2.7] that the terms 
in the above formula (for arbitrary type) correspond to chains of Weyl group elements giving 
rise to paths in the quantum Bruhat graph [2]; these chains are partitioned into segments 
corresponding to the parts of r, and the down steps (in Bruhat order) are measured by a 
statistic called level. In [311 we defined crystal operators on the above chains, both classical 
ones fi, i > 0, and the affine one /q. We called this construction the quantum alcove model, as it 
generalizes the alcove model in |321I33| . The main conjecture is that the new model uniformly 
describes tensor products of column shape KR crystals, for all untwisted affine types. The 
conjecture is proved in the same paper in types A and C, by showing that the bijections in 
[50] from the objects of the quantum alcove model to tensor products of the corresponding KN 
columns are affine crystal isomorphisms. With S. Naito and M. Shimozono, we are working on 
a uniform proof of the conjecture for all untwisted affine types. This will imply that the level 
statistic in the Ram- Yip formula expresses the corresponding energy function. In particular. 



we would obtain a generalization of the combinatorial formula (1.1) to arbitrary type, where 
the charge is replaced by the level. 

7.4. Demazure crystals and Macdonald polynomials. Ion yL2j showed that, in simply- 
laced types, one can identify a Macdonald polynomial at t = 0, namely Pfj,{x]q,0), with an 
affine Demazure character corresponding to a shift in the affine Weyl group. This raises the 



question of what happens in other types. The combinatorial formula (1.1) for P^{x;q,0) of 
type C, together with the decomposition of the corresponding tensor product of KR crystals 
with only the Demazure arrows (studied in Section l5| suggest that, in general, Pfi{x; q, 0) could 
be identified with a sum of Demazure characters. If so, one would want to determine them 
explicitly and, in particular, to know if all correspond to shifts in the affine Weyl group. 

As Ion's result and the Ram- Yip formula were given for non-symmetric Macdonald poly- 
nomials as well, we can also ask about the generalization of Ion's result to arbitrary types in 



the non-symmetric case. There are indications that the combinatorial formula (1.1), which was 



conjectured to have a version for arbitrary type (see Section 7.3), should be replaced with one 
whose right-hand side is a sum over a subset of the corresponding tensor product of KR crystals. 
If so, it would be desirable to describe this subset in an explicit way. See also [121 Section 9.2]. 
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